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Abstract 



If neutrinos produced in cosmological sources oscillate, neutrino telescopes 
can have a chance to detect signals from r-neutrinos. We present an estimate of 
^ ; the 'double bang' event rate produced by z/ r 's in km 3 scale detectors. 

m _ 

O ! 1. Introduction 

CO 

■ Up v T oscillations should lead to the proportion v e :v ii :v T = 1:1:1 for 

neutrinos produced in cosmological sources that reach the Earth, though the fla- 
vor ratio at production is expected to be 1:2:0 (if all muons decay). Passing 
through the Earth, u/s generate r-leptons via charged current (CC) interac- 
tions. Being short-lived particles, r's decay in flight producing z/ T 's and (in ~35% 
cases) also z/^'s or z/ e 's. Hence, neutrinos of all flavors undergo a regeneration 
process. The Earth is transparent for VHE v T 's and, moreover, z/ r 's are able 
^ ■ to recover partially also v^ e fluxes for which the Earth is opaque at VHE. i/ r 's 

can be detected via r's produced in CC interactions. At energies E v < 10 6 GeV 
r range is short (R T < 30 m) and showers at CC interaction and r decay ver- 
tices cannot be reconstructed separately in underwater /ice neutrino telescopes 
(UNT). Hence the 2 hadronic or hadronic+electromagnetic showers cannot be 
distinguished by a point-like hadronic+electromagnetic shower, produced by a v e . 
For E v ~ 10 6_; " 7 GeV the 'double bang' (DB) signature (r track between 2 showers 
at v T interaction point and r decay vertex [7] ) can be a clear and background- free 
evidence of v T detection. At E v > 3 • 10 7 GeV r-lepton range becomes larger than 
1 km and r's cannot be distinguished from muons in km 3 scale UNT. Using the 
muon and r propagation code MUM [13], we estimate that even though r's and 
muons have different energy loss properties, a lower energy muon can produce 
through stochastic processes showers with similar features than a r-lepton. 

Estimates of DB rates in km 3 UNT were done in [1] for down-going u/s, 
while events from the lower hemisphere were not considered. Calculations of up- 
going u T propagation were done in [3]. It was concluded that, though VHE is T 's 
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Fig. 1. Spectra of astrophysical z/s sampled according to the spectrum in [12] trans- 
formed after propagation through the Earth at 4 different nadir angles. Left column: 
no oscillations, v e :Vn'-V T = 1:2:0, right column: oscillations with maximal mixing, 
v e :v„:v T = 1:1:1. Thin solid lines: total incoming flux of astrophysical z/s; thick 
solid (right panels), dashed and dash-dotted lines (overlapping in right panels): 
out-coming v T + z/ T , + u^, v e + v e after propagation through the Earth, respec- 
tively. v^ e spectra include secondary neutrinos produced by u T 's. In the right part 
of each panel v atm spectra (conventional -l-prompt [10]) after propagation through 
the Earth are shown, upper thick markers: ZAj + lower thin markers: u e + T> e . 

are not absorbed by the Earth, their energies decrease; hence the amount of r 
events in the DB energy range is low. In this work we present results of a MC 
simulation of astrophysical v propagation through the Earth using 2 spectra [9,12] 
not considered in [1,3]. The rate of DB events in km 3 scale UNT is estimated. 

2. Method and Results 

We simulate both DIS CC and NC interactions for all kind of neutrinos 
[8] in the energy range 10 3 GeV < E v < 10 9 GeV using the CTEQ3_DIS structure 
functions. MUM [13] was used for r propagation including the newest corrections 
for photo- nuclear interaction [2]. We used TAUOLA package [6] to generate the 
r decays. Comparison of results for power-law spectra obtained by our algorithm 
with ones published in [4] showed a reasonable agreement. v^ e tracking is stopped 
after a CC interaction, r's produced in v T CC interactions are propagated up to 
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Fig. 2. The same as in Fig. 1 but astrophysical neutrinos are generated according to 
Mannheim, Protheroe and Rachen upper limit on diffuse neutrino flux [9]. 

their decay vertex. Incident fluxes on Earth of z/'s are assumed equal to v fluxes 
for all flavors. The background of atmospheric (ttK) and prompt (RQPM) i^e's 
is generated using spectra in [10]. Simulation results are presented in Fig. 1. We 
used the v flux predicted by Protheroe [12] for an external photon optically thick 
proton blazar model. In Fig. 2 results are presented or an upper bound (not 
for a model) on diffuse v spectrum for optically thick AGN sources (thick solid 
line in Fig. 5 in [9]). v T fluxes exceed z/^g ones remarkably up to nadir angles 
9 ~ 60° since z/ r 's are not absorbed by the Earth in contrast to z/^e's. But their 
spectra are shifted to lower energies due to energy degradation in regeneration 
process. For all 8, the outgoing flux of astrophysical neutrinos exceeds v a tm flux at 
E u > (1 -=-3) • 10 5 GeV. Secondary z/ M>e 's which are produced by r decays contribute 
to total outgoing u^ e spectra with fractions 0.57, 0.18, 0.06, 0.02 (spectrum [12]) 
and 0.18, 0.06, 0.02,' 0.01 (spectrum [9]) for nadir angles 6 =6°, 60°, 80°, and 85°, 
respectively. Analysis of ratio between muon rates and shower rates that may 
be an indirect signature of v T appearance was presented in [3]. In this work we 
analyze direct v T detection through DB signature. The rate of totally contained 
DB events is given by: 

/0(l) rco 
/ V eff (E l/T ,d)I(E UT ,e)a cc (E Ur )dE UT d(co S 9), (1) 

where Na is the Avogadro number, p is medium density (we use p = 1 g cm -3 
which is close to sea water/ice density), I(E Ut ,6) is the differential v T flux. The 
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Earth shadowing effect is accounted for the lower hemisphere (upper limit of 
the integral 0, while cos 6* = 1 is used to compute the number of events for the 
whole sphere); V e ff = S p (9) (L — R t (E Ut )), with S p projected area for tracks 
generated isotropically in azimuth at the fixed 9 directions on a parallelepiped 
(IceCube-like lxlxl km 3 [5] and NEMO-like 1.4x1.4x0.6 km 3 [11]), L geometrical 
distance between entry and exit point, R T r-lepton range, a cc total CC v cross 
section. E min = 2 • 10 6 GeV corresponds to r-lepton range R™ tn = 70 m. DB 
can be considered as a background-free r signature [7] since there should be 
no atmospheric fi which could produce 2 showers with comparable amount of 
photons. Nevertheless, this reasonable assumption needs to be verified through a 
full simulation. Table 1 shows DB rates for lower, upper and both hemispheres. 
Values for upper hemisphere are 3-r-6 times lower compared to [1] since more 
optimistic predictions for diffuse neutrino fluxes are used there. 



Table 1. Number of totally contained DB events in km 3 detector per year. 



Spectrum 


IceCube-like (JV_ 27r / N 2n / N i7T ) 


NEMO-like (AL 2?r / N 2n / jV 47r ) 


[12] 


0.7 / 1.6 / 2.3 


1.0 / 2.1 / 3.1 


[9] 


1.0 / 2.3 / 3.3 


1.4 / 3.1 / 4.5 



3. Conclusions 



Based on optimistic models [9,12] for AGN diffuse neutrino fluxes we found 
that in km 3 UNT one can expect a marginally observable rate of contained 'double 
bang' r events. 
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